Typeset with jpsj2.cls <ver.l.2> 



Letter 



o 
o 

(N 
>^ 



o 
o 

u 

Oh' 



X3 

o 

> 
oo 

(N 
O 

o 

-I— > 



I 

o 
o 



X 
S3 



Quasi-Classical Calculation of the Mixed-State Thermal Conductivity 
in s-Wave and d-Wave Superconductors 
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Department of Physics, University of Montenegro, Podgorica 81000, Montenegro 

To see how superconducting gap structures affect the longitudinal component of mixed-state 
thermal conductivity Kxx{B), the magnetic-field dependences of Kxx(B) in s-wave and d-wave 
superconductors are investigated. Calculations are performed on the basis of the quasi-classical 
theory of superconductivity by fully taking account of the spatial variation of the normal 
Green's function, neglected in previous works, by the Brandt-Pesch-Tewordt approximation. 
On the basis of our result, we discuss the possibility of Kxx{B) measurement as a method of 
probing the gap structure. 

KEYWORDS: quasi-classical theory, thermal conductivity, mixed state, s-wave superconductor, d-wave su- 
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1. Introduction 

In the last two decades, several new superconduc- 
tors showing unconventional pairing properties have been 
found among heavy-fermion, oxide, and organic materi- 
als."'^ These discoveries have motivated studies on identi- 
fying the gap symmetry, as well as clarifying the mech- 
anism of the superconductivity. At present, it has been 
well established that bulk measurements of mixed-state 
thermodynamic quantities can provide important in- 
formation on the gap structure. The principle be- 
hind this technique is that the low-energy quasiparti- 
cles around vortex cores, relevant to the low-temperature 
thermodynamics, sensitively reflect the structure of the 
superconducting gap. On the basis of the numerical so- 
lution of quasi-classical Eilenberger equations, we have 
demonstrated that the gap structure can be inferred from 
the low-temperature field dependence*'^ and field-angle 
dependence^*^"^^ of the specific heat and magnetization. 

Thermal-transport measurement in a superconduct- 
ing mixed state is another option for such a spectro- 
scopic method based on the bulk property of a sam- 
ple. At very low temperatures, where the transport 
property is mainly determined by the elastic impurity 
scattering, the thermal conductivity can also provide 
important information on the pairing state. Roughly 
speaking, the longitudinal component of mixed-state 
thermal conductivity Kxx{B) at zero temperature is 
given by k,,{B) ~ (i;27Vp(0)T^'^(0)), where iVp(O) is 
the momentum-resolved density of states at the Fermi 
surface, v^,. is the x-component of the Fermi velocity, 
Tp{0) is the momentum-dependent transport lifetime, 
and (• • • ) denotes the average over the Fermi surface 
(see eq. (30) below). As Volovik-^ pointed out, Np{0) 
and Tp{0) in a nodal superconductor are modified by 
the field-induced Dopplcr shift on the delocalized quasi- 
particles, while in an s-wave superconductor, such an ef- 
fect is negligibly small. Hence, we expect that the be- 
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havior of Kxx{B) in a nodal superconductor is different 
from that in an s-wave superconductor. Recently, a series 
of mixed-state thermal-transport experiments have been 
extensively performed by Matsuda's group to determine 
the gap structures in several newly found superconduc- 
tors.^^ 

The low-temperature mixed-state thermal transport is 
mainly determined by the following three mechanisms: 
impurity scattering, giving the Drude thermal conduc- 
tivity in the normal state; Andreev scattering by vor- 
tex cores; and the Doppler shift due to the supercur- 
rents."'^^ So far, there have been several approaches to the 
calculation of the mixed-state thermal conductivity in 
a moderately clean superconductor. I'^^^s Among them, 
two approximations have been frequently used. One is 
the Doppler shift approximation.^^' ^'^ This method ne- 
glects the Andreev scattering by vortex cores, and thus is 
valid only near the lower critical field Hd ■ The other ap- 



proximation 



is based on the Brandt-Pesch-Tewordt 



(BPT) approximation.-^^' The latter is superior to the 
former in that it takes account of both the Doppler shift 
and the vortex-core scattering. In the BPT approxima- 
tion, however, the spatial variation of the normal Green's 
function is neglected, which is strongly related to the lo- 
cal density of states detected by the scanning tunneling 
microscopy experiment. Hence, this approximation is 
valid only near the upper critical field i?c2- At present, 
there is no theoretical work valid in the wide field range 
from Hci to Hc2- In view of analyzing the experimen- 
tal data for several newly found superconductors, it is of 
importance to calculate the mixed-state thermal conduc- 
tivity beyond these two approximations. 

It is well known that the quasi-classical theory of su- 
perconductivity^^' ^'^ can accurately describe the mixed- 
state properties in a wide field range from Hd to Hc2- 
The advantage of this framework is that, if applied to 
the calculation of mixed-state thermal conductivity, it 
can take account of the spatial variation of the normal 
Green's function neglected in the BPT approximation, 
as well as the Andreev vortex scattering neglected in the 
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Doppler shift approximation. Of course, we need to solve 
a set of transport-like equations to accomplish the com- 
putation. 

In this work, we adopt the quasi-classical theory of su- 
perconductivity, and develop a method of calculating the 
mixed-state thermal conductivity Kxx{B). Then, we ap- 
ply our analysis to two-dimensional s-wave and d-wave 
superconductors, and calculate the magnetic-field depen- 
dences of Kxx{B) to clarify the effect of the gap structure 
on the mixed-state thermal transport. We also study the 
effect of the in-plane Fermi surface anisotropy. On the 
basis of our result, we discuss the possibility of using 
i^xx{B) as a method of probing the gap structure. 

The paper is organized as follows. In §2, we develop 
a method of calculating the mixed-state thermal con- 
ductivity based on the quasi-classical theory of super- 
conductivity. In §3, we present our numerical results for 
the mixed-state thermal conductivity obtained by solv- 
ing the quasi-classical equations. We discuss our result in 
§4, and conclude in §5. We use the unit h = c = = 1 
throughout this paper. 

2. Formulation 

2.1 Linear-response- equation 

Our approach to mixed-state thermal transport is 
based on the Kubo formula for a homogeneous tempera- 
ture gradient.^'* Then, the mixed-state thermal conduc- 
tivity Kxx{B) is given by 



Kxx{B) = -^Im 



(1) 



where Q^.^. is obtained from the (Matsubara) heat- 
current correlation function 



Qxxi^l^n 



dre" 



^{Tr[Jx{ri,T)Jxir2,0)])^ 



,(2) 



after the analytical continuation Q^^{u}) = Qxx{^<jJm 
u) -\- iO+). Here, u)m = "^nTm are bosonic Matsubara fre- 
quencies, (• • •)cq denotes the statistical average, and V 
is the volume of the sample. The heat-current operator 
J in the real-time representation is given by 



■1 e 



\A{2)) 



+ {Vr,-i\e\A{l))dt, VU^)V'(^) (3) 

where ^jJ is the electron-field operator and 1 — 
(ri, ti), 2 = (r2, t2)- To evaluate Qxx, we follow the pro- 
cedure suggested by Klimesch and Pcsch."^^ They showed 
that the heat-current correlation function can be con- 
structed in a manner similar to obtaining the density 
correlation hmction D under charge disturbance. On the 
imaginary (Matsubara) axis, D is given by 



dre"- 



x( - (T^[n(ri,r)n(r2,0)]) 



(4) 



where n — X^spinV'^V' is the electron density. The pro- 
cedure by Klimesch and Pesch is justified as long as we 
consider a moderately clean superconductor ^o/l <C 1 
(coherence length and mean free path /), where the 
impurity vertex corrections for the heat current can be 
safely neglected. When we consider a dirtier supercon- 
ductor, we have to include the vertex corrections to sat- 
isfy the conservation of the energy density. A full 
quasi-classical treatment for the charge current that in- 
cludes the impurity vertex corrections can be found in 
ref. 39. Note that the present result can be also derived 
from the Keldysh technique.^" 

Let us briefly review the procedure used by Klimesch 
and Pesch^^ to find the connection between the heat- 
current correlation function Qxx and the density corre- 
lation function D. Consider an external scalar potential 
$ = —\e\(j) that couples to the electron density n. Next, 
we define the quasi-classical Green's function as 



9 = 





' G 








-o) 



(5) 



90 



where G, F, and F"^ are the normal and anomalous 
Green's functions in the Gor'kov formalism, and is 
the single-particle energy measured from the Fermi en- 
ergy. We divide this Green's function into two parts, 
9 = 90 + Sg, where 

50 -i/o 
.i/o -50 

is the static part, and 

91 fi 

.-/2 52 

is the linear- response part. Then, using the Kubo formula 
for the charge disturbance, we obtain an expression for 
D within the quasi-classical accuracy, 

(8) 



(6) 



(7) 



£)(iw™) =i27r7VFT> (51 



where gi = .91/$, and N-p is the density of states in 
the normal state. Here we neglected the contribution far 
from the Fermi surface. On the other hand, by a di- 
rect evaluation of the density correlation function in the 
Gor'kov formalism, we obtain, within the quasi-classical 
accuracy, 



42 



D{iujm) = i2-kNfT 



^{GoGo-FoF^}), (9) 



in which we introduce a shorthand notation Go Go ~ 
FqF^ = Go(i£« + iwm)Go(ie„) - Fo{ien + iuJrn)F^{ien)- 
Go, -foi and Fq are the static components of the Green's 
function in the Gor'kov formalism. Note that, as men- 
tioned before, we neglect the impurity vertex corrections 
in this work. We also neglect the vortex motion, and con- 
sider only the quasiparticle contribution. Then, by com- 
paring the above two equations, we obtain the following 
relation. 



g,= f^(GoGo-FoF^ 

J ITT 



(10) 
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Next, we apply the same calculation scheme to the 
thermal conductivity Kxx{B). By a direct evaluation of 
the heat-current correlation function in the Gor'kov for- 
malism, we obtain within the quasi-classical accuracy. 



2.2 Numerical procedure 

To obtain an input quantity for the linear-response 
equation, eq. (13), we have to solve the static equa- 
tion (15). To do so, we adopt the Riccati parameteri- 
zation,45,46 



X J ^(GoGo-Fo4))(ll 



fo = 



2a 



Jo 



2b 



1-ab 



.90 



(18) 



From the last two equations, we can express Q^x in terms 
of 6g as 



51 



(12) 



The quantity Sg is determined by the linear-response 

equation:'*^''*"' 

V ■ V5g + En+TsSg - EnSgrs + i \e\v ■ Ats - A,i 
+ ^{9o{i£n+))Sg - Sg^{go{ien+)) 



go{ien+)i^ - i$ffo(ien), 



(13) 



where £„+ = £„ -|- w^, [A, = AB — BA, and 
T3 = (J.^*!). In the above equation, we assumed Born 
scattering for convenience, and unitarity scattering will 
be discussed later. Eq. (13) is supplemented by the nor- 
malization condition, 



go{ien+)Sg + Sggo{ien) = 0. 



(14) 



The static part of the quasi-classical Green's function 
go satisfies the following static (or Eilenberger) equation. 



e„T^ + {go)/2T+i\e\v A-iA,go =0, (15) 



with the normalization condition = 1- Here, A(p, r) = 

i-A*(r) ^'0''^)' ^^'^ {" ') denotes the Fermi surface av- 
erage. 

Once 5g is obtained, Kxx{B) can be calculated from 

cqs. (1) and (12). The Matsubara-sum in cq. (12) can be 
transformed into a real frequency integral using the well- 
known formula T^^^ /(i^n) = / ^th(^)/(£). Then, 
we obtain 

, 2 



l + ah' 1 + ab' 1 + ab' 

by introducing two functions a and b. Then, the static 
equation is transformed into the following Riccati-type 

equations: 

vVa = -2(e„ + i|e|v- A)a-f Ap- A;a2,(19) 

vVb = 2{en + i\e\vA)b-A*p + Apb^, (20) 

where e„ = e„ + {g)/2T, Ap = Ap + (/o)/2t, and r is 
the quasiparticle mean free time. The pair potential is 
expressed as Ap(r) = 'WpA{r) with the order parameter 
A(r) and the pairing function Wp. For a d-wave (s-wave) 
superconductor, we use Wp = V2{Px ~ Py)/PF (li'p = 1)) 
with pp being the Fermi momentum. 

The numerical procedure used to solve eqs. (19) and 
(20) is described in ref. 47. Calculations are performed 
for a two-dimensional superconductor with a hexagonal 
vortex lattice state. As mentioned before, we consider 
a moderately clean superconductor ^0/^ = 0-1, where 
^0 = v-p/2-kTco and I = vpr, with vp being the Fermi 
velocity. We initially assume a two-dimensional isotropic 
Fermi surface, and later we discuss the effect of the Fermi 
surface anisotropy. We also assume that the supercon- 
ductor has a large Ginzburg-Landau parameter kgl ^ 1 
(most of the newly found superconductors satisfy this 
condition), so that the magnetic field B is assumed to be 
constant. 

First, we determine the order parameter A(r) self- 
consistently by the following gap equation. 



Kxx{B) = ANpT 



de 



(2t) /,.2r :~-Rl 



/o 2rcosh2(^) 

(16) 

where gi^is) = giiisn — > £ + iO+), and [• • -Jsp denotes 

the spatial average. If we perform the low-temperature 
expansion of the above expression, we have 



= -Np{vl[-ig,%^). (17) 

T— »0 O 

Finally note that, as we will show later, the compact 
result^^ due to the BPT method corresponds to an ap- 
proximate expression of eq. (16). 



In ( A(r) ^2nTJ2( M) 



A(r) 



(21) 



where Tco is the transition temperature at zero field with- 
out any impurity. This procedure requires the anoma- 
lous Green's function /o, and it is obtained by solving 
eqs. (19) and (20) on the Matsubara axis. Note that for a 
d-wave superconductor, we set the anomalous self-energy 
(/o) equal to zero. This treatment is consistent with the 
neglect of the impurity vertex corrections for the heat 
current, since both quantities vanish in the limit of zero 
magnetic field. Of course, for an s-wave superconductor, 
we cannot neglect (/o), which satisfies the Anderson the- 
orem for pair breaking. 

In Fig. 1, we show real-space images of the obtained or- 
der parameter A(r) for an s-wave superconductor. Near 
Hc2 (Fig. 1(a)), vortices overlap with each other. How- 
ever, upon lowering the field H <^ Hc2 (Fig- 1(b)), the 
vortex core radius becomes smaller and smaller than 
the intervortex spacing. These behaviors can be ob- 
tained only after we calculate the order parameter self- 
consistently. After we determine the order parameter 
A(r), we solve eqs. (19) and (20) for the real frequency 
i£„ — > £ ± iO+ to obtain the retarded (advanced) Green's 



4 J. Phys. Soc. Jpn. 



Letter 



Author Name 



|A(r) I 




(a) 



B~Hc2 

(B/Hc2=0.765) 



|A(r) I 




(b) 



B«Hc2 

(B/Hc2=0.085) 



Fig. 1. Real-space images of the order parameter in an s-wave 
superconductor, obtained from our self-consistent calculation. 
r = {x,y) and A(t') are normalized by = {2|e|B)~^/^ and 
Ao = 1.764Tco, respectively. |A(t')| in a high field [(a)] and a 
low field [(b)] are shown. 



function. The impurity sclf-cncrgy is determined self- 
consistently in this static equation. '^'^ 

Next we solve the linear-response equation. To calcu- 
late Kxx{B) using cq. (16), we need to carry out the 
analytical continuation of eq. (13) to real frequencies. 
The off-diagonal components of 5g can be eliminated by 
the normalization condition, eq. (14). Then, the diago- 
nal components of the linear-response equation, eq. (13), 
yield 



(L + M) 



91 

92 



i$(<7«-5-4) 



Here, L and M are defined by 
L = 



V ■ V 





M = 



a + a 
a 





-V ■ V 

/? 

a + (3 



(22) 



(23) 



(24) 



where a{r) = ^"^0^+^}^° ^ Ap/^ +a;/„^ ^ 

^ ' aa^9o ' 90 -90 ' ^ ' 

- 9^), and Ap = Ap + ^{fi')_ 
To solve eq. (22), we first divide M into two parts. 



M = [M]sp -I- SM, and rewrite eq. (22) as 

91 \ _ ;^t„R „An / 1 



(i+[M]sp) 



92 



i<3>(5^-5^) 



1 



SM 



91 

92 



(25) 

Since gi and g2 are gauge-invariant and periodic, eq. (25) 
can be solved iteratively by a Fourier transform.^^ 

2.3 Relation to the BPT approximation 

For an analytical treatment, it is convenient to in- 
troduce the BPT approximation, which has been 
applied to the approximate calculation of the thermal 
conductivity. ^^'^^ Recently, this BPT approximation has 
been used frequently to interpret experimental data.^*^' 
To derive the compact result using the BPT approxi- 
mation, we first neglect the spatial variations of gi and 
g2 in eq. (22) and remove L. Next, we perform a spa- 
tial average on each matrix element of M as a(r) 

[Ap/t"+A;/o-^]3p .. . [Ap/t^+A;/«]3p 

'h^[{9a' ^ 9o)]sp- Note that the spatial average is per- 
formed separately in the denominator and the numera- 
tor. Then, gi and 52 can be obtained by a matrix inver- 
sion, and this yields 

de (^)^ 



8NfT 



I 

Jo 



:Re[g 



2T cosh^ (^) 
x{vlRe[g^UT^'{e)), (26) 
Re[Ap4« + A;/o«], 



Re[g^ 



^2l) 



|sp 



If we substitute the result of the BPT approximation for 
the static components of the quasi-classical Green's func- 
tion, then the above expression reproduces the compact 
result derived in ref. 25. 

3. Numerical Results 

First, we show our numerical results for the density of 
states. The density of states in a superconductor is given 

by 

N{e) 



Re[(.9o"(e))],p. 



(28) 



Of particular interest is the zero-energy density of states 
N{e = 0), i.e., the density of states at the Fermi en- 
ergy. This is because N{e = 0) is related to the low- 
temperature specific heat C through N{e = 0)/Nf = 
C/Cn\t^o, where Cn is the normal-state specific heat. 

In Fig. 2, we plot the zero-energy density of states 
in s-wave and d-wave superconductors {^o/l = 0.1) ver- 
sus magnetic field B. As in the impurity-free case,^ we 
can see the well-known difference in their field depen- 
dences even in the moderately clean superconductor. For 
an s-wave superconductor, we find a linear dependence 
of N{e = 0)/Np on B, which originates from the local- 
ized quasiparticles within the vortex cores, while for a 
d-wave superconductor, we find an approximate VB de- 
pendence due to the field-induced Doppler shift on the 
delocalized quasiparticles. A feature specific to the im- 
pure case is that there is a residual value of N{e = 0) in 
the zero-field {B — > 0) limit. 
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2D, T/Tco=0. 1,^0/1=0.1 




o l . . . . 

0.2 0.4 0.6 0.8 1 

B/Hc2 

Fig. 2. Magnetic- field dependence of zero-energy density of states 
N{e = 0) in s-wave (filled squares) and d-wave (filled circles) 
superconductors in the moderately clean case (^o/l = 0.1). Here, 
ATp is the density of states in the normal state. 



2D, T/Tco=0. 1,^0/1=0.1 




0.2 0.4 0.6 0.8 1 

B/Hc2 

Fig. 3. Magnetic-field dependences of the longitudinal thermal 
conductivity Kxx(B) in s-wave (filled circles) and d-wave (open 
squares) superconductors. KN,a;a: = ^ti^v'^Tt is the thermal con- 
ductivity in the normal state. Filled squares represent the zero- 
temperature value of Kxx{E) [see eq. (17)] in a d-wave supercon- 
ductor. 

Next, we show our numerical result for the thermal 
conductivity Kxx{B). Figure 3 shows the magnetic-field 
dependences of the calculated thermal conductivities in 
s-wave and d-wave superconductors. In the figure, we 
can observe a considerable difference in their field de- 
pendences: in an s-wave superconductor, Kxx{B) at low 
fields is exponentially small, while in a d-wave supercon- 
ductor, there is a finite amount of thermal transport even 
at low fields. This means that the field dependence of 
Kxx{B) can be used as a method of probing gap nodes. 



2D, ^o/l=0.1 




ol : . : . 1 

0.2 0.4 0.6 0.8 1 

B/Hc2 

Fig. 4. Magnetic-field dependence of zero-temperature thermal 
conductivity in a d-wave superconductor in the unitarity limit. 

However, we would like to emphasize here that unless 

wc discuss the zero-temperature limit, Kxx{B) in both 
superconductors have similar field dependences. This is 
evident when we look at the data in Fig. 3 at finite tem- 
peratures. Actually, it is difficult to distinguish Hxx{B) 
in a d-wave superconductor from that in an s-wave su- 
perconductor except for their residual values, although 
the data arc plotted in a moderately low temperature 
region {T/Tco = 0.1). Therefore, if we aim to use Kxx{B) 
to detect gap nodes, we have to extract its value at zero 
temperature. 

In some unconventional superconductors, thermal con- 
ductivity measurements in zero magnetic field sug- 
gest the necessity of including the resonant impurity 
scattering in the unitarity limit. ^^"^^ This effect can 
be taken into account by performing the replacement 
57 (.9o) ^(.9o)~"^ in eqs. (13) and (15). Figure 4 shows 
the magnetic-field dependence of zero-temperature ther- 
mal conductivity in a d-wave superconductor with the 
unitarity-limit scattering. The result suggests that the 
magnetic-field dependence of Kxx{B) in the unitarity 
limit is slightly different from that in the Born limit, 
particularly at low fields, while the main conclusion of 
the previous paragraph remains unchanged. 

Here, we comment on the difference between the 
present result and that of the BPT approximation. In 
the BPT approximation, the spatial average is performed 
incoherently as noted above using eq. (26). This proce- 
dure is expected to overestimate the vortex-core scatter- 
ing and reduce the thermal conductivity. Actually, if we 
compare Fig. 1 of ref. 25 with the present result, we can 
see that Kxx{B) obtained by the BPT approximation is 
slightly suppressed to a lower value. It has weak field de- 
pendences in the intermediate fields, and shows a steep 
increase near Hc2- In the present study, in contrast, we 
do not use such an artificial spatial averaging procedure, 
and we take account of the spatial dependence of Green's 
function in both static and linear-response equations. 
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2D, ^0/1=0.1 




0.2 0.4 0.6 
B/Hc2 



Fig. 5. Magnetic-field dependences of Kxx in an s-wave super- 
conductor witli an anisotropic Fermi surface. Tiie inset shows 
the corresponding data for N(e = 0)/Np. The definition of the 
anisotropy parameter b is given in eq. (29). 
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d-wave superconductors. 



A carcftil reader may think that the inchision of Fermi 
surface anisotropy may strongly affect the result. In- 
deed, we have realized that an in-plane Fermi surface 
anisotropy does affect the ficld-anglc-dcpcndcnt specific 
heat and magnetization oscillations.^^' To consider this 
possibility, we investigate the effect of in-plane Fermi 
surface anisotropy on the magnetic-field dependence of 
i^xx{B). Let us consider a model dispersion. 



zero-temperature expression of eq. (27), 



ep = ^(Px+Py)(l + &cos(4^)), 



(29) 



where (p denotes the azimuthal angle of Fermi momentum 
p, and the parameter h represents the degree of in-plane 
Fermi surface anisotropy with a fourfold symmetry. Here, 
we consider a case where the superconducting gap is an 
isotropic s-wave-type but the Fermi surface has in-plane 
anisotropy. 

Figure 5 shows the magnetic-field dependences of 
KxxiB) for b = 0.2. Note that the degree of anisotropy 
(6 = 0.2) is sufficiently strong to reverse the in-plane 
magnetization oscillation pattcrn'^^ in a d-wavc super- 
conductor. Concerning the field dependence of Kxx{B), 
there is a difl[erence between results for & = and 6 = 0.2. 
However, the difference is quantitative, whereas the dif- 
ference between s-wave and d-wave results (Fig. 3) is 
qualitative. This means that the superconducting gap 
structure influences the field dependence of Kxx{B) more 
strongly than the Fermi surface anisotropy. 

4. Discussion 

We discuss the origin of the difFerence in Kxx 

{B) be- 
tween a rf-wave superconductor and an s-wave supercon- 
ductor; a d-wave superconductor can transport a certain 
amount of heat even at low fields, whereas in an s-wave 
superconductor, the thermal transport is completely sup- 
pressed. For this purpose, it is convenient to employ the 



T ^ . ,--.v-x— L.u.=.-pv- 0)). (30) 

In Fig. 6, we show the magnetic-field dependence of 
{Tp{e = 0)) in both superconductors. By comparing 
this figure with Fig. 3, we can conclude that the main 
difference in Kxx{B) between s-wave and d-wave super- 
conductors results from that in the transport lifetime 

We can explain this difference in the following way. 
In an s-wave superconductor, the Andreev scattering 
rate given by the second term on the right-hand side of 
eq. (27) increases rapidly upon decreasing the field due 
to the fact that it contains the small factor Re[5Q-]sp ~ 
B/Hc2 in the denominator. In a d-wave superconductor, 
quasiparticles moving along the antinodal direction are 
dominated by the same principle. However, nodal quasi- 
particles do no_t undergo^ the Andreev scattering mecha- 
nism since Rc[Ap/Q^ + A* /(f ]sp is almost zero, and they 
can transport a certain amount of heat. It is this effect 
that yields the difference (B) between a nodal su- 

perconductor and a fully gapped superconductor. 

This interpretation is consistent with the following 
physical picture provided by Boaknin et al.^^ In an s- 
wave superconductor, localized quasiparticles within vor- 
tices, which are relevant to the low-temperature proper- 
ties, can give a linear-field density of states, but cannot 
contribute to the heat transport due to their localized 
nature. In a d-wave superconductor, on the other hand, 
Doppler shifted nodal quasiparticles can contribute to 
the heat transport, as well as to the well-known ^/B- 
dependent density of states. Thus, Kxx{B) can pick up 
only the contribution of extended quasiparticles. To some 
extent, this picture can be inferred from "extended zero- 
energy density of states A'^ext(5)" defined in ref. 9 (see 
Fig. 1 therein). Our study confirmed this physical picture 
by a concrete numerical calculation. 
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In this work, wc only consider a single-band supercon- 
ductor with a two-dimensional Fermi surface. Hence, we 
do not present a detailed comparison with experiments. 
Instead, we discuss the implications of the present results 
on experiments. For s-wave superconductors, the present 
result is able to explain the following experimental ob- 
servations. The magnetic-field dependence of Kxx{B) in 
Fig. 3 has a good correspondence to the experimental 
data for Nb (Fig. 3 of ref. 56) and VsSi (Fig. 3(c) of 
ref. 57), although our result cannot be applied directly to 
the Iow-kgl superconductor Nb. Recently, Kasahara et 
al. also observed a similar magnetic-field dependence of 
i^xx{B) in a pyrochlore superconductor KOs206.^* Con- 
cerning the d-wave calculation, our result is consistent 
with the thermal conductivity of Bi2Sr2CaCu208 at low 
fields. If we limit ourselves to the low-field region, our 
result is also consistent with the experimental data in 
other (perhaps nodal) superconductors,^^' in that 
these data at low fields have stronger field dependences 
than the s-wave superconductors discussed above. On 
approaching the upper critical field, however, a marked 
difference appears. Our result shows {Hc2 — S)^/^ be- 
havior near Hc2, which is consistent with the previous 
work,^'' while the experimental data show {Hc2 — B)^ 
dependences with r}>l. We leave this discrepancy to a 
future work. 

5. Conclusion 

In this paper, we have studied the magnetic-field 
dependence of the mixed-state thermal conductivity 
i^xx{B) in s-wave and d-wave superconductors based 
on the quasi-classical theory of superconductivity. The 
present work is valid in a wide field range from Hd to 
Hc2, and beyond the previous theoretical studies; the 
Doppler shift approximation is valid only near Hd , and 
the Brandt-Pesch-Tewordt approximation is valid only 
near Hc2- 

On the basis of our result, we clarified that there is 
indeed a clear difference in the low-temperature mixed- 
state thermal transport between s-wave and d-wave su- 
perconductors; in a d-wavc superconductor, there is a fi- 
nite amount of thermal transport even at low fields, while 
in an s-wave superconductor, the thermal transport is 
completely suppressed. This suggests that the magnetic- 
field dependence of Kxx{B) can be used to distinguish 
a nodal superconductor from a fully gapped supercon- 
ductor, on the condition that we extract a suitable zero 
temperature-limit. Also, we microscopically confirmed 
the physical picture provided by Boaknin et al.^^ on the 
mixed-state thermal transport, and demonstrated that 
the above-mentioned difference in Kxx{B) mainly results 
from the corresponding difference in the transport life- 
time. 

Acknov^rledgements 

One of the authors (H. A.) would like to thank H. 
Ebisawa, Y. Matsuda, T. Shibauchi, and Y. Kasahara 
for discussions, and M. Eschrig for important comments. 
He is also grateful to Y. Nisikawa and A. Furukawa for 
introducing him to the application of Kubo formalism to 
thermal transport. A part of the numerical calculations 



was carried out on Altix3700 BX2 at YITP in Kyoto 
University. 



9 
10 
11 

12 
14; 

18 
19 
20 

21 
22 

23 

24 
25 

26 

27; 

28 

29 

30 
31 

32 
33 

34 

35 

36 



37 

38 

39 

40 
41 



M. Sigrist and K. Ueda: Rev. Mod. Phys. 63 (1991) 239. 

C. CaroU, P. G. de Gennes, and J. Matricon: Phys. Lett. 9 

(1964) 307. 

G. E. Volovik: JETP Lett. 58 (1993) 469. 

H. Won and K. Maki: Europhys. Lett. 30 (1995) 421. 

N. B. Kopnin and G. E. Volovik: JETP Lett. 64 (1996) 690. 
S. H. Simon and P. A. Lcc: Pliys. Rev. Lett. 78 (1997) 1548. 
L Vekhtcr, P. J. Hirschfeld, J. P. Carbottc, and E. J. Nicol: 
Phys. Rev. B 59 (1999) R9023. 

M. Ichioka, A. Hasegawa, and K. Machida: Phys. Rev. B 59 
(1999) 184. 

N. NaJjai, P. Miranovic, M. Ichiolsa, and K. Machida: Phys. 
Rev. B 70 (2004) 100503(R). 

P. Miranovic, M. Ichioka, N. Nakai, and K. Machida: Phys. 
Rev. B 68 (2003) 052501. 

P. Miranovic, M. Ichioka, N. Nakai, and K. Machida: J. Phys.: 
Condens. Matter 17 (2005) 7971. 

H. Adaclii, P. Miranovic, M. Ichioka, and K. Machida: Phys. 
Rev. Lett. 94 (2005) 067007. 

H. Adachi, P. Miranovic, M. Ichioka, and K. Machida: J. Phys. 
Soc. Jpn. 75 (2006) 084716. 

For a review, see Y. Matsuda, K. Izawa, and I. Vekhter: J. 
Phys.: Condens. Matter 18 (2006) R705-R752. 
A. F. Andreev: Zh. Eksp. Teor. Fiz. 46 (1964) 1823 [Sov. Phys. 
JETP 19 (1964) 1228]. 

K. Maki and T. Tsuneto: Prog. Tlicor. Phys. 27 (1962) 228. 

K. Maki: Phys. Rev. 158 (1967) 397. 

R. M. Cleary: Phys. Rev. B 1 (1970) 169. 

A. Houghton and K. Maki: Pliys. Rev. B 4 (1971) 843. 

F. Yu, M. B. Salamon, A. J. Leggett, W. 0. Lee, and D. M. 
Ginsberg: Phys. Rev. Lett. 74 (1995) 5136. 

M. Franz: Phys. Rev. Lett. 82 (1999) 1760. 

C. Kubert and P. J. Hirschfeld: Phys. Rev. Lett. 80 (1998) 

4963. 

Yu. S. Barash and A. A. Svidzinsky: Phys. Rev. B 58 (1998) 
6476. 

H. Won and K. Maki: cond-mat/0004105. 

I. Vekhter and A. Houghton: Phys. Rev. Lett. 83 (1999) 4626 
(see Fig. 1 therein for comparison). 

M. Takigawa, M. Ichioka, and K. Machida: Eur. Phys. J. B 27 
(2002) 303. 

S. Dukan, T. P. Powell, and Z. Tesanovic: Phys. Rev. B 66 

(2002) 014517. 

A. C. Durst, A. Vishwanath, and P. A. Lee: Phys. Rev. Lett. 
90 (2003) 187002. 

U. Brandt, W. Pcsch, and L. Tewordt: Z. Phys. 201 (1967) 
209. 

W. Pesch: Z. Phys. B 21 (1975) 263. 

H. F. Hess, R. B. Robinson, R. C. Dynes, J. M. Valles, Jr., and 
J. V. Waszczak: Phys. Rev. Lett. 62 (1989) 214. 

G. Eilenberger: Z. Phys. 214 (1968) 195. 

A. I. Larkin and Yu. N. Ovchinnikov: Zh. Eksp. Teor. Fiz. 55 
(1968) 2262 [Sov. Phys.-JETP 28 (1969) 1200]. 
R. Kubo, M. Yokota, and S. NaJjajima: J. Phys. Soc. Jpn. 12 
(1957) 1203. 

R. Klimesch and W. Pesch: J. Low. Temp. Phys. 32 (1978) 
869. 

In the Meissner state (B = 0), a Ward identity ensures that 
the impurity vertex corrections vanish for even-parity super- 
conductors.^^' 

V. Ambegaokar and L. Tewordt: Phys. Rev. 134 (1964) A805. 
P. J. Hirschfeld, P. Wolfle, and D. Einzel: Phys. Rev. B 37 

(1988) 83. 

M. Eschrig, J. A. Sauls, and D. Rainer: Phys. Rev. B 60 (1999) 
10447. 

P. Miranovic: unpublished. 

G. M. Eliashberg: Zh. Eksp. Teor. Fiz. 61 (1971) 1254 [Sov. 
Phys. JETP 34 (1972) 668]. 



8 J. Phys. Soc. Jpn. 



Letter 



Author Name 



42) K. Scharnbcrg: J. Low. Temp. Phys. 6 (1972) 51. 

43) A. L Larkin and Yu. N. Ovchinnikov; J. Low. Temp. Phys. 10 
(1973) 407. 

44) Yu. N. Ovchinnikov: Zh. Eksp. Teor. Fiz. 66 (1974) 1100 [Sov. 
Phys. JETP 39 (1974) 538]. 

45) N. Schopohl and K. Maki: Phys. Rev. B 52 (1995) 490. 

46) N. Schopohl: cond-mat/9804064. 

47) P. Miraiiovic, M. Ichioka, and K. Machida: Phys. Rev. B 70 

(2004) 104510. 

48) P. W. Anderson: J. Phys. Chem. Sohds 11 (1959) 26. 

49) P. Miranovic and K. Machida: Phys. Rev. B 67 (2003) 092506. 

50) L. Tewordt and D. Fay: Phys. Rev. B 64 (2001) 024528. 

51) H. Kusunose, T. M. Rice, and M. Sigrist: Phys. Rev. B 66 
(2002) 214503. 

52) C. J. Pethick and D. Pines: Phys. Rev. Lett. 57 (1986) 118. 

53) P. J. Hirschfeld, D. Vollhardt, and P. Wolfle: Sohd State Com- 
mun. 59 (1986) 111. 

54) S. Schmitt-Rink, K. Miyake, and C. M. Varma: Phys. Rev. 



Lett. 57 (1986) 2575. 

55) E. Boaknin, R. W. Hill, C. Proust, C. Lupien, L. Taillefer, and 
P. C. Canfield: Phys. Rev. Lett. 87 (2001) 237001. 

56) J. Lowell and J. B. Sousa: J. Low. Temp. Phys. 3 (1970) 65. 

57) E. Boaknin, M. A. Tanatar, J. Paglione, D. Hawthorn, F. Ron- 
ning, R. W. Hill, M. Sutherland, L. Taillefer, .1. Sonier, S. M. 
Hayden, and J. W. Brill: Phys. Rev. Lett. 90 (2003) 117003. 

58) Y. Kasahara, Y. Shimono, T. Shibauchi, Y. Matsuda, S. 
Yonczawa, Y. Muraoka, and Z. Hiroi: Phys. Rev. Lett. 96 
(2006) 247004. 

59) H. Aubin, K. Behnia, S. Ooi, and T. Tamegai: Phys. Rev. Lett. 
82 (1999) 624. 

60) H. Suderow, J. P. Brison, A. Huxley, and J. Flouquet: Phys. 
Rev. Lett. 80 (1998) 165. 

61) K. Izawa, Y. Nakajima, J. Goryo, Y. Matsuda, S. Osaki, H. 
Sugawara, H. Sato, P. Thalmeier, and K. Malci: Phys. Rev. 
Lett. 90 (2003) 117001. 



